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A theoretical investigation of population dynamics in an A^-level system driven by pulses of arbi- 
trary shape is presented, which is based on Floquet theory. The pulse only has to have finite support 
in time, which is well fulfilled in experiments. Furthermore the dynamics must be approximately 
time-local, which restricts the applicability to negligible memory effects. The rotating-wave approx- 



qh' imation is not invoked. We discuss control of both the real and imaginary part or just the modulus 

Cn . °' ^"^ ^°"P°"^"' °' ''^ '"^^ wavefunction. A necessary criterion on the Floquet quasienergies 

^ , of pulses of arbitrary shape is obtained for population inversion in the case N = 2, which is a 



generalization of the vr-pulse criterion previously demonstrated by Holthaus et. al. We find that 
only the eigenvectors of the propagator contain enough information to give a sufficient criterion. 



^^ . By analysing the propagator in an A^-level system we estimate the number of real parameters of 

the pulse that have to be adjusted to control transitions between arbitrary initial and final states, 

l^ ' including superposition states. The propagator can be recast into block diagonal form, leading to 

j^ ■ a description by an effective 2-level-system composed of the initial and final state, for which our 

criterion is still valid, however additional conditions on the pulse are required. In all cases the 
number of conditions on the pulse parameters increases linearly with A^. 
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I. INTRODUCTION 

For a long time there have been attempts to control specific quantum transitions with 
high efficieny using pulses that might itself consist of a train of subpulses [l|. While a 
number of approaches rest on combining experiment and simulation using feedback 2|], in 
the present paper we study model systems in the most general setup to find simple criteria 
of control. There are many applications where such criteria would be desirable, as in charge 



excitation and transport in biological systems [3|, constructing efficient QBits for quantum 



a 



a 



computers |^, or spin injection and magnetic transport in spintronic devices 5|. 

Not even the two-level system (2LS) can be solved analytically for the general case. 
Traditional approximations are the rotating wave approximation (RWA) |6| and Floquet 
theory [7|. The former is justified only for small detunings between field- and transition- 
frequency while the latter in its original form is applicable for strictly period fields (where 
it is exact), although in adiabatic approximation it can be extended to sufficiently slowly 



a. 



varying pulse envelopes 

In recent years short pulses of only few optical cycles [9| have become available, which are 

not suitable for adiabatic approximations. For such pulses semiclassical strong field theory 

12| applies, which however is valid only in the limit of large quantum numbers. Therefore 

most theoretical studies have used numerical techniques, including the relatively new topic 



of dynamics in molecular systems under the infiuence of such strong pulses [lOj. Recently 
we contributed an approximate treatment of A^-level systems in the presence of such pulses 

Regarding full population inversion, Holthaus and Just have previously used Floquet 



theory to analyse population transfer in a pulse with adiabatic envelope within the RWA 13 
Il4| . In the present paper we show analytically that a generalization of that criterion provides 
a necessary condition for population inversion under arbitrary pulses of finite support under 
fairly general conditions, i.e. without invoking RWA (retaining however the semiclassical 
dipole approximation). An important additional observation is that a sufficient condition 
for population inversion is provided only by the Floquet eigenvectors. We generalize the 
discussion by investigating either full control of the wavefunction (populations and phases) 
or just the populations. The former case is important for interference experiments. We 
evaluate the number of conditions on real parameters of the pulse, noting that the nonlinear 



equations for the pulse parameters resulting from the Floquet equations might have no 
physical solution. 

We add a detail concerning notation. Consider the setup of arbitrary orthogonal initial 
and final vectors in an A^-level system. Whenever the square modulus of each component of 
the fina wavefunction (at the end of the pulse) is equal to the corresponding quantity of a 
given final target vector, we say that the corresponding pulse achieves (generalized) "pop- 
ulation inversion" (PI); the case where initial and final vector are both system eigenstates 
corresponds to the usual notion of a resonance. If the wavefunction itself is equal to the 
final vector, we use the term "phase selective population inversion" (PSPI). 

On the other hand in case initial and final vector are not orthogonal we deal with a 
general control scenario and use the terms "population control" (PC) and "phase selective 
population control" (PSPC). As this scenario goes somewhat beyond the topic of population 
inversion, we discuss it separately in appendix |Al 

II. THEORY 

In the following we consider the interaction of a quantum system, approximated by A^ 
levels, with a laser pulse in terms of the semiclassical dipole approximation, noting that the 
derivations do not depend on the details of the interaction. Setting h = 1 and choosing a 
basis Bo := {b^^^ . . . , b^^^}, the time-dependent Schrodinger equation becomes 

idtipj = '^Hjkipk , Hjk = ej6jk - fijkE(t) . (1) 

k 

Here Sjk stands for the Kronecker symbol. The ipj, j = 1,. . . ,N, are the time dependent 
expansion coefficients of the state described by the wavefunction ip, the ej are the eigenvalues, 
and the Hjk are elements of the Hamiltonian H. E{t) is the projection of the field onto the 
dipole operator with expectation values fijk- The diagonal elements fijj represent permanent 
dipole moments. Note that our derivation also applies to non-hermitean Hamiltonians, which 



are used in studies of open systems [l5| , as long as memory effects are negligible during the 
pulse. 

In order to use the Floquet method we consider the experimental pulse E{t) of duration 
T to represent one period of an auxilary periodic field Ep(t). By the assumed time-local 
behaviour populations at time T will be equal for both fields. 



A. Floquet Theory 

We define the propagator transporting the basis B(to) = Bq as B(t) := U(t, to)B(to) and 
rewrite the Schrodinger equation as 

zG(t)U(t,to)=0 , G{t):=dt + tli{t) , (2) 

using the contravariant operator G with its local representation G(t) along the one- 
dimensional (Id) submanifold of propagation in the unitary group U(N). After completing 
a closed loop in t we obtain B(T + to) = U(T + Iq, to)B(to). In mathematics U(T + t^, to) 
is called the monodromy matrix. Note that when transporting along the loop for a sec- 
ond time we start with B(T + to). Due to the periodicity of the auxilary field Ep{t) 
we have G(t) = G(T + t), the propagator is equal to the one of the first loop, and 
B(2T + to) = [U(T + to, to)]^B(to). We conclude that the set of all loops forms an abelian 



group 16|, which then has only Id irreducible representations [IZ]. In this way the following 



exponential homomorphism is obtained, 

U(nT + to, to) = exp(OTriT) , flj^ = ujj6jk ■ (3) 

The Uj are the well-known Floquet quasienergies, which are real whenever H = H^ and 
complex otherwise. 

B. Population inversion in two- level systems 

In order to find criteria for PI it suffices to analyse the general form of the propagator 
matrix after a full period. We set to = for the time when the external perturbation is 
switched on, so that U(0, 0) = I2, the 2d unit matrix (below, I„ will denote the n x n 
unit matrix). Furthermore we abbreviate V := U(T, 0). The propagator of an A^-level 
system belongs to the manifold U(N). For A^ = 2 we expect 4 parameters to determine the 



propagator [17|. The equations VV^ = VtV = I2 lead to 

e^^° cos{6) e*^" sm{5) 
-e'^-sm{6) e*T'*cos(5) 

with the additional relation 7a + 7d = 7b + 7c + 2J7r, j G Z. Defining the quantities 



V 



7a, 76, 7c, ^yd^S eR (4) 



la + ld . la -Id . lb- Ic ,-. 

X ■■= ^ , Ai := 7\ , A2 := — - — (5) 



and extracting a phase factor we finally arrive at 



V = exp(ix) 



-e-*^2 sin(5) e-^^i cos((5) 



(6) 



The matrix exp(— 'ix)V has unit determinant and hence is an element of SU(2). Only the 
parameter 6 determines the modulus of elements of V, whereas all other parameters just 
contribute phases. 

Using a relation for the Wronski determinant from the theory of ordinary differential 
equations [l8|, namely 

det U(t, 0) = det U(0, 0) exp \ -i [ ti B.{t')dt' I , (7) 



we can immediately identify 



X 



tr U{t')dt' 



(8) 



We proceed with the setup of PI investigated in [13| and take as the initial state ip^^^ = 
[1,0]""" (a superscript T denotes the transpose, not to be confused with the period T) with 
population only in level i, and use ip^^^ = [0, 1]^ with population only in level f. In order to 
find situations more readily realizable in experiment we concentrate on the case n = 1 [n the 
number of pulses) and assume that such a pulse exists (as suggested by numerical results, a 
rigorous proof of existence for the general case is yet to be given). From the structure of the 
propagator in eq. ([6]) and tp we obtain l-ip^^^T)]"^ = sm"^ 5 for the population in f, reflecting 
the fact that with regard to PI only 6 is relevant. PI then demands 

2k + 1 



c-PI 



-TT 



kez 



(9) 



We note that 2k + 1 can be interpreted as the number of Pis during the pulse. 

In order to find a necessary criterion based on the Floquet quasienergies we diagonalize 
V, that is, we solve the polynomial 



det(V-CI) = 



(10) 



resulting in 



a,2 



exp{ix) I cos{5) cos(Ai) ± [cos^(5) cos^(Ai) — 4exp( 



-«x 



1/2 



:ii) 



In case of PI we find that the two values Ci,2 are symmetric with respect to 0, namely 



Ci^2 = ±^exp(i-) . 
This gives the necessary condition 

{u^^ - ujI^)T = (2n + l)7r , neZ 



(12) 



(13) 



as a criterion for arbitrary pulses, analogous to the one based on an adiabatic Floquet 

n r 

approach in refs. [13|| and [M]. Clearly eq. ( IT3l) corresponds to an avoided crossing situation 
because without the off-diagonal contributions in V the eigenvalues would both be zero. 

With a view on applications of condition f lT^ . we note that in a pragmatic sense at first 
sight it does not seem to lead beyond the results obtained by integrating the Schrodinger 
equation numerically. However, efficient expansions exist for the Floquet determinant, which 
could be used to find the quasienergies up to the required approximation from a time- 
independent system of equations [7|. Moreover an analytic approach might answer the 
question of existence of a pulse effecting full PI. 



Returning to eq. (ITTl) . we note that due to the vanishing trace of V, Ai 



2k+l 



71 also 



leads to eq. (TT3|) . This shows that this criterion itself does not provide a sufficient condition 
indicating PI. In order to find a sufficient criterion we investigate the eigenvectors e^ and 
e^, which are determined from the equation 



det(V - Cal)e° = . 
Using e^ = A^i[l, ci]^, e^ = A^2[c2, 1]"^ we obtain the equations 

Cl(2) 



Ci[C2 



sin((5) 



exp(— iA2) — exp[i(Ai — A2)] cot((5) 



(14) 



(15) 



Inserting the corresponding eigenvalue we arrive at 



1 



ci(c2^) = 2^^pH(^2 -x)]^ cot(5[cosAi -2exp[2(Ai -x)]] 



+ - 



cot^ S cos^ Ai — 



4exp(-ix) 
sin^ 6 



1/2- 



(16) 



In case of PI we obtain for the rotation matrix the simple expression 



,Pi 



V2 



1 — exp(ia) 

exp(— za) 1 



, a := A2 - -(x + tt) ,. 



(17) 



which clearly is unitary and provides the sufficient criterion for PI. 

We observe that if we set 6 = 6^^ in eq. (Q, the value of Ai is irrelevant. This is in 
analogy to spherical coordinates because they are a subset of the Euler angles parametrizing 
rotations in 0(3), which is homomorphic to the group SU(2) determining the propagator if 
X is fixed. Therefore we may interpret Ai as an azimuth angle while 6 would represent the 
inclination. The rotation matrix e^^ is real whenever A2 = and x = (2^ + I)?!", A; G Z. In 
this case Ci,2 are real and 001^2 is 0, vr/T. 

In control theory the wavefunction at the end of the pulse is prescribed and the inverse 
problem is solved |l|. In order to discuss PSPC in this context, we take -ip^^^ = [a,bj^, 
a, 6 G C In case of PSPC we thus get 4 relations describing the real and imaginary parts of 
a and b, 

a = exp[i{x + Ai)] cos{6) , 6 = — exp[z(x — A2)] sm{5) , (18) 

of which only 3 are independent due to |ap + |6p = 1. Therefore PSPC determines all 3 
parameters Ai, A2 and 6, as expected, while x is given by eq. (j8]). On the other hand PSPI 
only requires to manipulate the 2 parameters 6 and A2 due to the fact that Ai is irrelevant 
in case cos((5) = 0. 

C. N levels 

First we analyse the geometrical properties of the propagator by using the Floquet rep- 
resentation. Assuming the quasienergies to be real of the form Uj = QjUq, with Uq := 2n/T, 
and restricting to the ffist Floquet zone, qj G [0, 1), we distinguish the cases of (i) rational 
Qj = rij/m (with n^, m G N relative prime), and (ii) at least one of the qj is irrational. In the 
ffist case, after m pulses we arrive back in the initial state. Due to the translation structure, 
in terms of real parameters the m vectors lie on a great circle of the 2N — 1 dimensional 
unit sphere in C^. For m even we have the additional inversion property that to each state 
ip, its "mirror" image —ip is also among the m vectors. In the second case of irrational Qj 
the (infinite) number of vectors densely fill the great circle. 

Control of the dynamics is determined by the equation ^^^^ = Vip^^K In case of PSPC 
this results in 2N equations for the real and imaginary parts of the components of ip^^\ Due 
to normalization of ip^^'' we therefore have to fix 2N — 1 real parameters of the pulse for 
PSPC. 



If only the populations \ijjj p are relevant and the phases do not matter (setup of PC), 
we simply have to use the correct number of equations which amounts to replacing 2N by 
N. Using again ||'?/'*-^''|P = 1, this leads to A^ — 1 conditions for population control. 

As is well known, time- reversal symmetry is present if H(t) = H''"(— t) holds, and 
^(f)t .— ■?/;t(2") is mapped into ■i/'*-'^^ := ^"^{0) via U(0,T) = U^(T, 0) by the equation 
idtip\—t) = ■(/'"l"(— t)H'''(— t). When the field strength of the pulse is nonzero, this symmetry 
is not present because usually H(t) ^ H''"(— t). However, due to the translation symmetry 
of the Hamiltonian in T in the Floquet problem we have time-reversal symmetry when con- 
sidering only times that are multiples of T, due to H(nT) = H(— nT) = H^(— nT) with 
n G N if the system is not open. This symmetry does however not provide extra information 
on the number of conditions for PSPC, since the resulting equation t/'*-'''^ = ^jj^^'^^V can be 
directly obtained from ip^^^ = V^p^^^ by hermitean conjugation. In contrast, if ^^^'^ and ^p^^'' 
are orthogonal, an additional property is obtained, which leads to additional information as 
shown in the next section. 

D. Eflfective 2LS for orthogonal initial and final vectors 

Next we consider the simplest case of the A^-level propagator, namely that all population 
resides in a single level prior to the pulse, tpf^ = 6ik, and also after the pulse, tpl = 
exp(i/3)(5ffe. The dynamics in the present situation projects out column i of the propagator, 
^f = V,. 

In case A^ = 2 and i = 1 PSPI then leads to 5 = 6^^ and A™ = -\ J^ ti B.{t')dt' - /3 + 
{2k + 1)tt, A; G Z, that is only 2 parameters are relevant due to orthogonal initial and final 
vectors. The PSPI propagator for N = 2 reads 



V(2),PSPi = exp(,^) 



e^^^'" 



(19) 



For A^ levels this generalizes to Vf^k = for each k ^ i. This follows from V^i = for Z 7^ f and 
orthonormality of the A^ column vectors of V, which results in = Yli=ii ^.fc^zl = ~^f,fc^f*i 
for all k ^ i. 

Orthogonal initial and final vectors correspond to the case m = 4, where the inversion 
property and \J(T)ijj^^^ = ^^^^ as well as \J {—T)ijj^^^ = —tjj'^^^ hold. Using U(— T) = V''", this 



case leads to the special condition —ip^^^''^ = ^(')1'V. By the above argument, in row i of V all 
entries turn out to be zero except \4,f = — exp(— z/3). Using the same reasoning as above we 
find that in column f only the i-th element can be nonzero. This results in 2N — 3 additional 
(real) conditions, because \4i = has already been invoked above to enforce PSPI. 

For convenience we permute indices to get i= 1 and f= 2 such that the effective 2LS 
shows up in the PSPI propagator in terms of the block-diagonal form 



Y(iv),pspi 



Y(2),pspi 

v(^-2) 



(20) 



where V^^ ^^ is the propagator of the orthogonal subspace. It should be noted that 
Y(7V),pspi g u(2) © U(iV - 2) C U(A^) only at the end of the pulse. The propagator is 
not block-diagonal while the pulse is on. Relating to the Floquet quasienergies and eigen- 
vectors, we note that the 2x2 block is equal to V*^^)'^^^\ and hence eq. f lT3|) indicates 
possible PSPI for orthogonal initial and final vectors independent of A^. 

In case of an initial superposition state we can use a rotation to transform ^*^'\ tp^^^ to 
the first and second basis vector, respectively. In the rotated frame the propagator becomes 
Vr = RVR'''. Defining ip-R := Hipy the above derivation of the number of conditions extends 
to general initial states whenever —ip^ = '4'r^r holds. However, this clearly follows from 
—ip(^) = \'^ip(^) under rotation. We note that V^^^^ will be block-diagonal but V^^^^ will 
not in general assume such a form. 

We summarize that in the case of PSPI for arbitrary orthogonal wavefunctions an effective 
2LS is obtained by fixing AN — 4 real coefficients of the pulse. Thus the number of conditions 
may be larger than the A^^ — 1 parameters that determine the propagator. In this case we 
have to manipulate the /«// propagator to obtain an effective 2LS. Investigating A^ = 2 as an 
example we have 4N — 4 = 4, while the propagator is determined by only 3 parameters. The 
reason lies in the fact that when we demand Vu = and set the phase in V12, the second 
column of V is already determined. 

We note in addition that control is independent of the initial vector and only depends on 
the pulse, due to unitarity of the propagator and linearity of the Schrodinger equation. The 
propagator V therefore performs the same unitary transformation for every initial vector. 

If only populations matter we again have to replace 2N by A^. The same arguments as 
mentioned above in this section lead toA^ — 1-|-A^ — 2 = 2A^ — 3 conditions in case of a PI. 



Whereas non-orthogonal initial and final vectors correspond to a control scenario rather 
than to population inversion, we can use similar arguments to estimate the number of 
parameters necessary for this extended control. In appendix [X] we show that for this setup 
the number of conditions increases to 8A^ — 13 (PSPC) and 8A^ — 15 (PC) for generating 
a 2LS, due to the absense of any symmetry. In appendix [B] we show furthermore that 
depending on the circumstances the number of conditions required may be smaller than 
derived above. The reason lies on the one hand in the upper bound of A^^ — 1 parameters for 
SU(A^), which becomes essential for small enough A^. On the other hand exceptional values 
of parameters can make other parameters redundant. An example is the irrelevance of Ai 
when cos(5) = in case of A^ = 2. 

E. The adiabatic limit 

Previously Holthaus [l3|, [ij] proposed the following relation among the instantaneous 
Floquet quasienergies in case of PI by a pulse with adiabatic envelope, 

[ dt[n['"''''''\t) - nl'"''''''\t)] = {2n + 1)7: , neZ , (21) 

Jo 

by replacing exact energy eigenstates (no field) by integrals over time-dependent quasiener- 
gies in phase- factors, using Floquet theory for continuous waves [l3| . In a later publication, 
RWA on resonance (zero detuning) is compared with Floquet theory and the relation be- 



tween quasienergies calculated for this special case is used to conjecture that eq. (l2Til remains 
valid also for RWA with nonzero detuning [ij]. In the following we show how eq. fl2T|) can 



be rigorously deduced from eq. flT3|) whenever the envelope can be treated adiabatically. 

Assume a pulse of duration T, represented by the product of a periodic function and a 
slowly varying envelope. Note that in general the periods can even have different length, dp, 
with p labelling the period, if attached e.g. at zeros to give a continuous field. Due to slow 
variation of the envelope we make the assumption that the envelope is constant within each 
period p. We can represent the propagator V as a product of contributions over all periods, 

p 

v = Jjepc-^^^'^'^'^et (22) 

p=i 

using standard Floquet theory within each of the P periods. 0*^^^ denotes the matrix of 
Floquet quasienergies for period p and Gp transforms to the diagonal frame. Making the 

10 



well-known additional (adiabatic) approximation that 0^0p_i ~ In for every p, we arrive 
at V ^ Gpexp (-i J2p=i O^^^cip) e[. This has to be compared to ee~'^^Q\ with Q and 9 
corresponding to the whole pulse (time interval [0,T]), within our approach of periodically 
repeating the pulse. Note 6^6p_i ~ Itv applied to every p means that the envelope was 
actually replaced by a rectangle of duration T. Therefore we get O = 0p, for every p, and 
after finally using that G^Op_i ~ In implies that 0*^^^ — O*^^"^-' is very small, we arrive at 

1 ^ 1 /-T 

rj = - ^ n^P^dp ^ - / rftf2(^^'^^)(t) . (23) 

p=i -^0 

Together with eq. ( IT3l) this leads immediately to eq. ( 12T|) . We note that in case of a pulse 

with slowly varying envelope and all dp equal (i.e., a well defined carrier frequency), eq. ( ETj) 

yields extra information on the transient dynamics. The cost, however, are the restrictions 

imposed by the adiabatic approximation. 

III. SUMMARY 

In the present paper we investigate control of the full final wavefunction (phase-selective 
population control, PSPC) or just final population control (PC) within the semiclassical 
dipole approximation. The analysis is valid for any pulse having finite support (of length 
T) in time if the equation of motion is local in time (no or only adiabatic memory effects). 

For orthogonal initial and final vector, population inversion (PI) in the 2-level system is 
governed by a single parameter (here denoted 6). In terms of Floquet quantities we obtain 
the criterion {uf—UijT = {2n+l)TT, n G Z on the difference of Floquet quasienergies from the 
2d subblock of the propagator corresponding to the initial and final vector. This criterion is 
however only necessary, as it can be fulfilled by varying another parameter of the propagator 
not leading to PI. A sufficient criterion is provided by the shape of the Floquet eigenvector 
matrix. In contrast, for phase-selective population inversion (PSPI) 2 parameters of the 
propagator have to be set appropriately. Applying the adiabatic approximation we show 
that our criterion specializes to the expression presented by Holthaus et. al. [l3|, Il4| . 

For A^-level systems we find that for PSPC at most 2N real parameters of the propagator 
have to be set, while N conditions are sufficient for PC. The propagator at t = nT, n G N, 
can always be made block-diagonal with one subblock of dimension 2x2, that is, initial 
and final level form a 2d subspace, decoupled from the (A^ — 2)-dimensional subspace of the 

11 



remaining levels. In order to create such an effective 2LS in a given basis with preselected 
initial and final vectors, additional conditions must be fulfilled; for PSPI at most 4N — 3 
real parameters must be fixed while for PI 2N — 2 parameters are needed and moreover 
the criterion on the Floquet quasienergies derived for N = 2 holds true unaltered for any 
N > 2. In case of PSPC for non-orthogonal initial and final vectors we need to set 8A^ — 12 
real parameters to obtain the block-diagonal form with 2x2 subblock, while it is 8A^ — 14 
for PC. The numbers of conditions derived are in general upper bounds due to possible 
dependencies between parameters of SU(A^). 

Appendix A: Effective 2LS for non-orthogonal initial and final vectors 

We now treat the case that ■?/''''-' and ^^^^ are two non- orthogonal vectors of length 1, for 
which we derive the number of conditions necessary to obtain an effective 2LS. Remembering 
that the Floquet quasienergies can be written as Uj = qjUo, Uq := 27r/T and Qj = rij/m G 
[0, 1), we have in the present case m 7^ 4 or at least one of the Qj will even be irrational. The 
dynamics does not involve only 2 but an arbitrary number of all A^ levels. However, as all 

points of the trajectory are located on the great circle, it should be possible to transform the 

(2) 
propagator to block diagonal form, with a 2 x 2 block Vj^ corresponding to the subspace 

spanned by the initial and final vectors in a rotated frame. We first construct a basis 

consisting of ■?/'*■'•*, a second vector orthogonal to ■ip'^^^ and coplanar with ^^^^ and ^^^\ and 

N — 2 vectors orthogonal to both il)^^' and i})^^'. A rotation R relates the standard basis 

to this new one and the propagator becomes Vr = RVR^ Such a transformation always 

exists, and in the transformed basis at most the first and second levels are populated in the 

initial and final state. 

For control it remains to manipulate Vj^ . In contrast to orthogonal ip'^^'^ and ■i/''-^^ this 
requires setting all (complex) elements of Vr.j^, {j < 3} A {A; > 2} or {/c < 3} A {j > 
2}, to zero. Adding in case of PSPC the 3 parameters of Vj^ , we have to fix a total of 
2 X 2 X 2(A^ — 2) -|-3 = 8A^ — 13 real parameters to obtain an effective 2LS. The larger number 
of parameters is required due to the lack of any symmetry. 

In case of PC only one parameter, 6 of V^ , is relevant. However we also have to set 
the 8A^ — 16 appropriate elements of Vr to zero, and thus we still need to impose 8A^ — 15 
conditions. 



12 



Appendix B: Restrictions on the number of conditions 

As noted in subsection IIIDl the number of conditions derived for control and for creating 
an effective 2LS may exceed the total number of parameters of the propagator. Noting 
that the multiplicative phase x i^i gQ- dS]) is predetermined by tr(H), we require N"^ — 1 
parameters specifying the propagator. The number of conditions derived in subsection IIIDl 
and appendix lAl is relevant only if it is sma//er than A^^ — 1, which in case of PSPC leads to 
A^ > 2, and for an effective 2LS to A^ > 3 (PSPI) and A^ > 6 (PSPC). 



^loquet matrix (which 
7|. This adds one real 



In order to solve the homogeneous system of equations from the 
in general is infinite), the Floquet determinant has to be set to zero 
condition in case of a hermitean Hamiltonian. 

Assuming A^ > 1 this leaves in total 2A^ real parameters for PSPC, and min(4A^ — 3, A^^) 
or min(8A^ — 12, A^^) real parameters of the pulse to obtain an effective 2LS for orthogonal 
and nonorthogonal vectors, respectively. The total number of conditions are A^ for only PC, 
and 2A^ — 2 or min(8A^ — 14, A^^) to obtain an effective 2LS for PC with orthogonal and 
nonorthogonal vectors, respectively. For better overview we present results for A^ < 10 in 
table E 

There is a further possible reduction of the number of parameters for necessary for control. 
For an illustration we come back to A^ = 2 and note that for orthogonal initial and final 
vector the parameter Ai was irrelevant in case cos((5) = 0. This situation of pure states 
(only one level of ip^^^ip'^^'' populated in the appropriate basis) restricts to a subset of SU(2). 
For A^ > 2 we expect such situations to occur as well, such that our numbers of conditions 
become upper bounds in general. Although a general parametrization of SU(A^) via a faithful 
matrix representation of its algebra with minimal dimension has been given |19J , more work 
is necessary to find detailed corrections to the number of conditions. 
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N 


Control 


Effec 

orthogonal ■i/'^'^^W 


tive 2LS 

non-orthogonal ^('),^(^) 


PSPC 


PC 


PSPI 


PI 


PSPC 


PC 


2 


4 


2 


3 


2 


4 


2 


3 


6 


3 


9 


4 


9 


9 


4 


8 


4 


13 


6 


16 


16 


5 


10 


5 


17 


8 


25 


25 


6 


12 


6 


21 


10 


36 


34 


7 


14 


7 


25 


12 


44 


42 


8 


16 


8 


29 


14 


52 


50 


9 


18 


9 


33 


16 


60 


58 


10 


20 


10 


37 


18 


68 


66 



TABLE I: Number of conditions as a function of the number of levels N for control and for 
generating an effective 2-level system for orthogonal and non-orthogonal initial ip^^^' and final tjj^ ' 
vector, respectively. Note due to depedencies of parameters of SU{N) the necessary conditions 
may be less than given here for N > 2. 
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